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Abstract 

Quantum liquids are characterized by the distinctive properties such as the low temperature 
behavior of heat capacity and the spectrum of low-energy quasiparticle excitations. In particular, 
at low temperature, Fermi liquids exhibit the zero sound, predicted by L. D. Landau in 1957 and 
subsequently observed in liquid He-3. In this paper, we ask a question whether such a characteristic 
behavior is present in theories with holographically dual description. We consider a class of gauge 
theories with fundamental matter fields whose holographic dual in the appropriate limit is given 
in terms of the Dirac-Born-Infeld action in AdS p +i space. An example of such a system is the 
J\f = 4 SU(N C ) supersymmetric Yang-Mills theory with Nf massless M = 2 hypermultiplets at 
strong coupling, finite baryon number density, and low temperature. We find that these systems 
exhibit a zero sound mode despite having a non-Fermi liquid type behavior of the specific heat. 
These properties suggest that holography identifies a new type of quantum liquids. 

PACS numbers: ll.10.Wx, 11.25.Tq, 61.20.Gy, 67.10.-j 
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Introduction. — Gauge/gravity duality [H, Q, 0| has become a useful tool for investigating 
strongly coupled field theories. In the class of models where this tool can be applied, the 
strong coupling limit of the field theory is mapped into the weak-coupling, classical limit of 
a gravity theory, which can be studied either analytically or with minimal computer power. 
For example, a cousin of QCD — the M = 4 supersymmetric Yang-Mills (SYM) theory — 
has been studied using this method. Such studies have pointed to a universal value of the 
viscosity/entropy density ratio in a wide class of strongly coupled theories (for a review, see 
0]). Somewhat surprisingly, the viscosity /entropy density ratio of the quark-gluon plasma 
created at the Relativistic Heavy Ion Collider seems to be close to this value, indicating that 
gauge/gravity duality may be useful for studies of QCD. 

Here we would like to see what the gauge/gravity duality has to say about strongly coupled 
quantum liquids. By quantum liquids we mean translationally invariant systems at zero (or 
low) temperature and at finite density. Given the important role that quantum liquids play 
in physics, it is natural to ask whether the newly developed technique of gauge/gravity 
duality can give us any insights into their behavior. 

The cornerstones of our understanding of quantum liquids are two phenomenological 
theories. These are Landau's Fermi liquid theory , \^ and the theory of quantum 

Bose liquids 0, 0] • These two theories describe two distinct behaviors of a quantum liquid at 
low momenta and temperatures. In a Bose liquid, the only low-energy elementary excitation 
is the superfluid phonon with a linear dispersion. This leads to a T 3 behavior of the specific 
heat at low temperatures. The Fermi liquid has a richer spectrum of elementary excitations, 
consisting of fermionic quasiparticles and a bosonic branch, which contains, in particular, 
the zero sound. The fermions dominate the specific heat, which scales as T at low T. 

In this paper, we found, through the gauge/gravity duality, a new type of quantum liquid. 
The quantum liquid we consider has a T 6 behavior (~ T 2p in p spatial dimensions) of the 
specific heat at low temperature. Despite the non-Fermi liquid behavior of the specific heat, 
the system supports a sound mode at zero temperature, which we will call "zero sound." 
The mode is almost identical to the zero sound in Fermi liquids: not only the real part of 
its dispersion curve is linear in momentum (u — v q) , but the imaginary part has the same 
q 2 dependence predicted by Landau a long time ago for quantum attenuation of the zero 
sound. The difference is that in our case the zero sound velocity coincides with the first 
sound velocity, while in the case of a Fermi liquid the two velocities are not equal to each 
other. 

In Fermi liquids, zero sound is a collective excitation involving fermions near the Fermi 
surface. It was predicted by Landau [6] and experimentally observed in liquid Helium-3. 
In a weakly interacting Fermi gas, the zero sound velocity is close to the Fermi velocity 
v F . The first sound, which is a hydrodynamic mode that exists at finite temperatures and 
wavelengths larger than the mean free path, has velocity v F /^/3 at weak coupling. The 
quantum attenuation (i.e. damping at zero temperature) of zero sound was first considered 
by Landau, who showed that the imaginary part of the zero sound energy scales as q 2 where 
q is the momentum. The experimental situation with the measurement of the zero sound 
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quantum attenuation is summarized in [K 

A specific example considered in this paper is the M = 4 SU(iV c ) supersymmetric Yang- 
Mills (SYM) theory with Nf massless M = 2 hyper mult iplet fields. This theory has been 
suggested as a model which approximates QCD better than the theory without fundamental 
quarks. A string-theoretic description of this system is given by a low-energy limit of the 
D3/D7 brane configuration. The theory has been studied at finite temperature and density 
using the gauge/gravity duality 11, 12, 13|, 0, [l5, 16, 17, 18, 19]. Nevertheless, two striking 



aspects of this and similar systems (characterized by the Dirac-Born-Infeld action in Anti-de- 
Sitter space) — the unusual behavior of the low-temperature specific heat and the existence 
of the zero sound — have so far eluded attention. Their description is the main purpose and 
the main result of the paper. 

Preliminaries. — The Af = 4 SYM theory contains fields in the adjoint representation of 
the gauge group only. Fields in the fundamental representation can be introduced by using 



the following construction 20j. In type IIB string theory, one considers a system of N c 
D3-branes and Nf D7-branes aligned in flat ten-dimensional space as 

Xo X\ Xi X3 X4 X5 Xq £7 X$ Xg 

D3 x x x x (1) 
D7 xxxxxxxx 

In the limit of large number of colors (N c ^> 1) and large 't Hooft coupling (g\ M N c ^> 1), the 
D3-branes are replaced by the near-horizon AdSs x S 5 geometry [Ij, while the Nf D7-branes 
can be treated as probes embedded into this geometry as long as Nf/N c <C 1, i.e. as long as 
their backreaction on the geometry can be neglected 
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The standard form of the near-horizon D3 brane metric is 

ds 2 = ^ n^dafdx" + ^ (dr 2 + r 2 d{l 2 5 ) , (2) 

where r]^ = diag(— 1, 1, 1, 1), R is the curvature radius of the AdS 5 (we shall set R — 1 in 
the following). The horizon is located at r = 0. 

In this paper, we focus on adding Nf massless M = 2 hypermultiplets to M = 4 SYM, 
keeping the theory at zero temperature. In the dual gravity this is described by a zero 
temperature "horizon-crossing" D7-brane embedding in which the distance between the D7 
branes and the horizon in the x§ — Xg direction vanishes (l3j ]. 

The action for the D7-branes is the Dirac-Born-Infeld (DBI) action 



Sdbi = -N f T D7 J d 8 £ v /-det(^ a6 + 27ra , F a6 



), (3) 

where T^ 7 is the D7-brane tension, £ a are worldvolume coordinates, g a b is the induced 
worldvolume metric and is the worldvolume U(l) gauge field. At the boundary, the 
gauge field couples to the U(1)b "flavor" current J M , where U(1)b is the "baryon number" 
subgroup of the global symmetry group U(Nf) possessed by the M = 2 hypermultiplet fields. 
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(The exact form of the U(1)b current operator is given in [ill].) Considering finite "baryon" 
density (J°) 7^ in the boundary gauge theory corresponds to turning on a non-trivial 
background worldvolume gauge field A (r) in the bulk hj. The DBI action becomes 



S 



DBI 



-Af Vz drr 



(4) 



where the factor 2na.' is absorbed into A , V 3 is the spatial volume of the boundary gauge 
theory, and the prefactor Af = \NfN c / (2n) A is determined by the gauge/gravity duality 
dictionary 
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The construction above is specific to the D3-D7 system, but we can be more general and 
consider Dg probe branes whose worldvolume include an AdS p+ 2 factor. For probe branes 
corresponding to massless flavors the embedding is independent of the internal directions. 



One example with p 
action reads 



2 would be the defect D5 on AdS 4 x S 2 in AdS^ x S E 
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S 



DBI 



-Afq V p 



dr r l 



1 — A' 2 
1 ^0 > 



The DBI 



(5) 



where the normalization now includes the tension of the Nf Dg-branes 
the embedding problem is given by 



131. The solution to 



A' 



d 



Vr 2 P + d 2 



(6) 



where d = {2na'N q )~ l p is proportional to the baryon number density p [13j. In all subse- 
quent formulas, the results for the D3-D7 case are trivially recovered by setting p = 3 and 
using the relation a'^ 1 = VX 

Low-temperature limit of the specific heat. — One interesting hint to the nature of the 
phase of matter described by the probe branes with finite chemical potential is the behavior 
of the specific heat at low temperature. To extract this information one first needs to 
generalize the setup to a background that contains a black hole, that is the AdS p+ 2 part of 
the metric is modified to 



p+1 \ 

^— ) dt 2 + dx 2 



+ 1 



'_H 



(7) 



Here r H is the horizon radius of the black hole related to the temperature of the black hole 
by r H = AttT/ (p + 1) (with R — 1). The action describing the system at finite density and 
finite temperature is almost identical to the one in the zero temperature case 



S 



DBI 



-M q V v 



drr p \l 



A' 2 



(8) 



the only difference being that the integration starts at r H rather than at 0. All powers of 
the redshift factor cancel between g u and g rr . The solution for A' Q is still given by Eq. fl^D. 



1 Throughout this paper, we work in a gauge with A r = 0. 
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The on-shell value of the total action 2 directly gives us minus the thermodynamic poten- 
tial Q in the grand canonical ensemble. Substituting the solution ([6]) for A' Q into the action, 
we can write Q = fi a( j + fifun, where f2 ad ~ T p+1 is the contribution of the adjoint degrees 
of freedom (for the D3/D7 system O ad = -7r 2 iV 2 T 4 /8 is the free energy of M = 4 SYM at 
strong coupling), and 

where A is the ultraviolet cutoff. In Eq. (Q, the local counterterm action built from the 
metric h^ u induced on the slice r = A by the ambient metric (J7D has been added in the spirit 



of the holographic renormalization [22[|. In the grand canonical ensemble, the potential Q 
is a function of T and fi, where \i is the baryon number chemical potential related to the 
density and temperature via the condition 

/■oo 

ji= / drA' . (10) 

Jr H 

The integrals in Eqs. flU]) and (TTDT) can be expressed in terms of the Gauss hypergeometric 
function: 

KV n r 2 J +1 „ (I . 1 1 i»\ . A^Kr? +1 



fifun - fi °"V/+^ 2Fl V2' 1 + V 2 + T P ] -f) + 2(p + l) ' (U) 

(11 1 \ 

2> ii+ v -*)■ (i2) 

where and /^o are the zero-temperature values, 

fi = adp, tt = - , (13) 

(p + l)a p 

and a(p) = V - ±\ T (l + ±\ /Y (\). We note that the last term in Eq. (fTTD . 

AfqVprp' 1 /2(p + 1) = O c .t., is independent of the matter density and has the same tem- 
perature dependence as the free energy of adjoint fields f2 ad . We shall focus on the density- 
dependent part of the thermodynamic potential AQ = Q {un — Q c t . Equations (iTTj) and (|T2|) 
determine AQ as a function of temperature and chemical potential. At low temperature, 
both equations can be treated as series expansions in T/[i <C 1. The baryon number density 
is proportional to d, 

'—vrfc-"'*- (14) 

One then computes the entropy density s(fi,T) in the grand canonical ensemble 

1 / dAQ(T,fi) 



•v^vA—ir*)^ (15) 



2 The total action is the sum of the DBI action ([8]) describing fundamental degrees of freedom and the bulk 
gravitational action dual to the adjoint sector of the system. 
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Using Eq. (fT2|) , we find the entropy density as a function of temperature and charge density 
s(T,d) = s +Af q [—-) — l + 0[Td-A , (16) 



p+lj 2d 

where so = 47ip/[(p + l)(27ra')] is the entropy at zero temperature. This entropy is related 
to the quark thermal mass (free energy), which is negative and proportional to T. Finally, 
the specific heat (heat capacity per unit volume) c v at constant volume and density is 
determined by 

At low temperature (T <C /x ) the density-dependent part of the specific heat 3 is proportional 
to T 2v : 

( 4tt \ 2p+1 T 2 p r / i\l 

c v =Kp\ 7 — l + 0[Td- P ) . (18) 



K p+lj d 

This has to be contrasted with a gas of free bosons whose low temperature specific heat is 
proportional to T p (a sphere of volume T v of occupied states in momentum space, each with 
energy T) or a gas of fermions, whose low temperature specific heat scales as T for any p 
(a shell of thickness T of occupied states above the Fermi surface contributing an energy T 
each). The behavior of the specific heat in Eq. ( fl~8l) is suggestive of a new type of quantum 
liquid. 

Zero sound. — The zero sound mode would manifest itself as a pole of the zero-temperature 
retarded flavor current density correlator 0, S, @]- In the dual gravity language, the pole 
arises as the quasinormal frequency of the background geometry 23, [2J, |25|. Generically, 
the quasinormal spectrum is determined by fluctuations of all background fields including 
the metric. However, in the particular case we are dealing with, it is sufficient to consider 
fluctuations of the DBI U(l) field in the gravitational background ([2]) with the non-trivial 
background component Aq. Moreover, since the dual quantum field theory is isotropic, we 
can choose the fluctuations to depend on time, radial coordinate and one of the spatial 
coordinates (e.g. x p ) only 

A^r) -> A M (r) + a tl (r,x , x p ) . (19) 

Substituting f|T9l) into the DBI action ([5]) and expanding to second order in fluctuations, 
we find that the quadratic part of the resulting action is given by the sum of the actions 
describing longitudinal (ao, a p ) and transverse fluctuations. The action for the longitudinal 
fluctuations is (we use the a r = gauge) 

S (2) = K j^ xdr J — (20) 

2 J 1(1 -^2)3/2 rVT^ yi^Vj 



The density- independent part of the specific heat is proportional to T p . 
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Introducing the Fourier components 



dui dq 



(21) 



we find the equations of motion for the fluctuations 



d 
dr 

d 
dr 



r p a' n 



:i-^o 2 ) 3/2 



r p-A 
1 — A n 



r p a' 



r-P-4 



+ /l- A 12 ^ qa ° + U2(X ^ 



0. 



(22) 



(23) 



There is also a constraint arising as a consequence of the residual gauge invariance of the 
components a and a p 

ua' + (l-A' 2 )qa' p = 0. (24) 
Introducing a new radial coordinate z — l/r, Eqs. (I2"2"l [231 |2"4"|) can be written as 

(25a) 
(25b) 
(25c) 



8 Z [f z 2 - p a' ) - f z 2 ~ p (uqa p + q 2 a ) = , 
3 Z (/ z 2 - p a' p ) + / z 2 - p [uqa + u\) = . 
f 2 uja' + qa = , 



where f(z) = y/1 + d? z 2p . Following the approach of 25], we use Eq. (125]) to derive the 
equation for the gauge- invariant variable E = uo a p + q qq: 



E" + 



f (Sq 2 - u 2 f 2 ) p - 2 
f(q 2 -uj 2 P) 



,, 2 f 2 a 2 

J 



(26) 



The horizon, z = oo, is an irregular singular point of the differential equation (1261) . The 
solution in the vicinity of z = oo is given by E(z) ~ e ±tuJZ /z. The incoming wave boundary 
condition at the horizon 24j singles out one of the exponents 



E(z) = C 



where C is a constant. For uz <C 1 we have 



(l + 0(l/z) 



E(z) = ^ + iuC. 



(27) 



(28) 



On the other hand, for wzC 1 and qz -C 1 with uj/q fixed, Eq. ( 1261) reduces to the equation 



E" + 



f (3q 2 - u 2 f 2 ) p - 2 



E' = 



/ (g 2 - a; 2 / 2 ) « 
whose solution is given in terms of the Gauss hypergeometric function 



E(z) = C x + C^" 1 



_ _ 2Fl U' 2 ~ V 2 ~ V " d ^ 



_Pf(z) Pip- 1 ) 



(29) 



(30) 
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For z — > oo we find 

E(z) -»• d + C 2 (- + &) + 0(l/z 2 ) , (31) 
where the coefficients a and 6 are given by 

tf-P^-r(j-j)rfe) 
a ~T' 6 ^r(l) ' (32) 

Matching to the expansion ( 1281 . we find the coefficients Ci and C2: 

d = fiw - - J C , C 2 = - . (33) 
\ ay a 

The lowest quasinormal frequency is found by imposing the Dirichlet condition at the bound- 



ary, E(0) = [251 ] . This condition gives the equation Ci = or, equivalently, 



d»r (3 - 5) r (4 



'"=i>-pj V6 • (34) 



Solving Eq. (13"4"|) for small u; and 0, we find the dispersion relation for the zero sound 

g ^ (|) g 2 

^ ^r(|-i)r(i) 



W = ±4= - : ^ + 0(g 3 ) . (35) 



Using the expression for the chemical potential at zero temperature from Eq. ( [TBI , the zero 
sound dispersion relation can be written as 



, ±J--ii- + Ouh- (36) 



What is the nature of this excitation? First, one can exclude the possibility that it is a 
superfluid phonon. Indeed, our background does not break the particle number symmetry, 
hence the ground state is not a superfluid. Furthermore, superfluid phonon width has a 
low- momentum behavior different from q 2 , namely q 5 in 3 spatial dimensions [7| and q p+2 
in p spatial dimensions (provided that phonon decay is kinematically allowed). The q 2 
behavior of the imaginary part is characteristic of zero sound quantum attenuation, thus we 
call this mode the zero sound. Yet in other respects (such as the specific heat temperature 
dependence) the system does not show Fermi-liquid behavior. It is notable that the zero 
sound velocity in our system coincides with the velocity of the finite-temperature first sound, 
while in a weakly-coupled Fermi liquid it is ^Jp times larger than the first-sound velocity. 

Conclusion. — In this paper we have considered a general theory described by a DBI 
action in AdS space. We found that by turning on a chemical potential one arrives to a new 
type of quantum liquid. The specific heat c v has an unusual non-Fermi liquid T 2p behavior 
(T 6 in 3+1 dimensions and T 4 in 2+1 dimensions). The low energy spectrum contains a 
gapless mode with a dispersion relation similar to the zero sound in Fermi liquids. One can 
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speculate that the mode observed here is what the Fermi-liquid zero sound becomes when 
the interaction is infinitely strong. In this connection, we note that in a simple model of the 
Fermi liquid, the velocity of the zero and first sounds approach each other in the limit where 
the interaction strength (parameterized by the Fermi-liquid parameter F ) is infinite [jj. 

The systems described here are strongly coupled, as they have gravity duals. It would 
be interesting to investigate the properties of the ground state and the zero sound in the 
weak-coupling regime of the M = 4 SYM theory with M = 2 matter hypermultiplets. We 
leave this problem for future work. 
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